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ABSTRACT 


The experimental set-up for laboratory study of spacecraft control/structural 
interaction has been designed. Design specifications have been derived, and all 
the actuators and sensors have been selected except the end-point displacement 
sensing of the arm. The mainbody and the flexible arm have been fabricated to 
meet design criteria. 

The equations of motion for the experimental model have been derived 
and natural frequencies determined. The natural frequencies of the fle.xible arm 
have been determined experimentally and compared with analytical predictions 
obtained by using the GIFTS finite element analysis program. The experimental 
and analytical results are in good agreement except the first mode. 
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I. INTRODUCTION 


A. THE CONTROL OF FLEXIBLE STRUCTURES 

Flexible structures have become an integral part of 
modern spacecraft design for a variety of reasons. 

Deployable structures extending several meters into space 
used for solar arrays are required for todays high powered 
satellites. Radar reflectors of large diameter are used for 
communication satellites. Meanwhile, as the systems expand 
in size, weight constraints continue to plague the satellite 
designer. Hence, lighter and flexible structures are often 
used. This results in lower structural natural frequencies. 

The large antenna reflectors require higher attitude 
pointing accuracy. In order to provide higher attitude 
pointing accuracy, the control bandwidth has to be 
increased. Because of the decrease in the structural 
natural frequencies and increase in control bandwidth, there 
is a high possibility of the structural frequency falling 
into the control bandwidth, resulting in control/structure 
interaction. 

For some space applications, the flexible structure will 
require active control. As a result, current and future 
spacecraft study of control/structure interaction has become 
a challenging problem. 
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B. BACKGROUND. 


The control of rigid body dynamic systems has been a 
subject of study for many years. Controlling three-axis 
stabilized satellites has been successfully achieved for 
several years. The 1990's, however, bring a new challenge 
to the control problem. Future NASA, military, and 
commercial space missions will involve advanced space 
systems which have higher power reguirements, greater 
reguired pointing accuracies, and faster slewing through 
larger angular motions. Additionally, permanent space 
stations, like Freedom, and complex space-platforms, like 
the Hubble Telescope, will reguire the use of large, 
flexible manipulators for maintenance, on-orbit 
construction, etc. 

Research into the demanding problems of the control of 
flexible structures has been undertaken for the last two 
decades by major space-system engineering firms, the 
government laboratories, and universities. The effort has 
been, however, limited to analytical studies. Recently, 
experimental work has started at some universities. 

Texas A&M University has constructed a hub-appendage 
configuration to perform large angle maneuvers with 
vibration suppression for a flexible space vehicle (Junkins, 
1989, pp. 1-4). The central hub pivots on a ceramic bearing 
in the horizontal plane, and four identical cantile'-’ered 
flexible appendages with endpoint masses protrude from this 
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base. A reaction wheel is driven by a DC brushless motor to 
provide a system torque. Sensors include an angle encoder, 
strain gauges and a motor tachometer. Control laws were 
successfully written to control the system; however, 
nonlinear friction/stiction occurred as a result of the 
ceramic bearing, and nonlinear bending moments were observed 
as a result of the flexibility which hampered accurate 
modeling of the system response. 

Old Dominion University has conducted an experiment to 
investigate the slewing of flexible structures while 
simultaneously suppressing vibrational motion during the 
maneuver (Yang, 1989, pp.l). The experimental setup 
consists of a trolley on which is mounted a long, flexible 
beam. The trolley is driven through a driver pulley and a 
cable transmission system. The flexible beam is rotated in 
a horizontal plane by the beam motor. Sensors include 
strain gages, angular potentiometers and a tachometer for 
the beam motor. Again, nonlinear effects were discovered in 
large bending deflections, and friction of the cables and 
trolley, but again they posed no problem in controlling the 
system. 

Stanford University has set up an experiment most similar 
to the one designed here, and mutual cooperation 
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between Stanford University and Naval Postgraduate School 
has proved invaluable in the early stages of design. 

A two link manipulator is used by Stanford University 
consisting of two flexible beams pinned at the shoulder, and 
joined at the elbow by a limited angle torquer motor 
(Oakley,1988, pp. 1-4). Rotary variable differential 
transformers measure joint angles, and a CCD television 
camera tracks endpoint position. The control equations were 
derived using an assumed modes method, and the experimental 
results concurred with the simulation. Again, despite 
nonlinear effects, accurate modelling was achievable. 

At the Naval Postgraduate School, analytical/ 
experimental studies have been started this year on the 
attitude control of flexible spacecraft. 

C. EXPERIMENTAL DESIGN 

The design of the experimental set-up requires a 
compromise between simulation realism and the practicalities 
of the hardware. A design is needed that will be an 
accurate model for understanding spacecraft control 
problems. For this reason, the experiment was constrained 
to rotate, not translate, about its principal axis. Pitch 
motion is the only motion being investigated and controlled. 
In space, a momentum wheel can be used to control pitch axis 
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motion since the pitch axis is the only axis that remains 
inertially fixed for Earth-oriented satellites. 

There are seven phases of development in the experimentc1 
design. This phased approach to building the experimental 
laboratory allows interim goals to be identified and 
realized, 

1. Phase I 

In the initial phase, the design will consist of a 
mainbody, an L-beam flexible arm, and a reflector (Figure 
1), The pitch motion wil^ be controlled by a motor driven 
momentum wheel. There will be an angle and angul'r rate 
sensor on the mainbody. The flexible arm end point position 
will be measured. Beam mode-shapes will be determined bv 
strain gages and accelerometers. 

The mainbody will be controlled with information from 
the ma'nbody only. The arm sensors will be utilized for 
performance measurements and mode shape analysis. The 
reflector will not be controlled during this phase. 

2. Phase II 

In phase II, tne central body will be controlled with 
information from the mainbody rate and angular position 
sensors and the endpoint position sensor. 

3. Phase III 

A limited-angle stepper motor will be added in phase 
III at the reflector-arm joint to control the orientation of 
the reflector with respect to a fixed point on the mainbody. 
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The angular position of the reflector with respect to the 
end point will be measured. The central body will be 
controlled, as above. The reflector will also be controlled 
in both open loop, closed loop with feedback from sensors on 
the reflector, and closed loop with feedback from all 
sensors. 

4. Phase IV 

Thrusters will be added to the mainbody to desaturate 
the momentum wheel. 

5. Phase V 

The strain gages and accelerometers will be fed into 
the control loop for controlling the reflector. 

6. Phase VI 

Liquid tanks will be added to the system. 

7. Phase VII 

A limited angle stepper motor and angle position 
sensor will be added at the L-joint and active control of 
the flexible arm will be achieved. Eventually, robotics 
experiments can be explored by fixing the mainbody. 


7 








D. OBJECTIVES 


This thesis is concerned with the initial design of the 
experiment through Phase T. There are three main 
obj ectives. 

1. Experimental Set-up 

The preliminary design of the experimental set-up 
includes the identification of system requirements, the 
resulting specifications, and the selection of sensors, 
actuators, and computer systems to make the experiment 
operational through Phase I. 

2. Equations of Motion 

The equations of motion for the experimental system 
must be written to determine the mass and stiffness 
matrices. 

3. System Identification 

The system characteristics (the natural frequencies 
and mode shapes), should be determined by computer aided 
finite eleinen* analysis and verified by experimental modal 
analysis. 
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II. EXPERIMENTAL SETUP 


The initial design of the experimental set-up consisted 
of determining the required specifications, researching the 
available equipment, and selecting the appropriate product. 
Whenever possible, the set-up was designed as a scaled-down 
version of existing space-platforms. The specifications 
were also often based upon real-world considerations, such 
as acceptable sizes and weights for a laboratory set-up. 

A. GRANITE TABLE 

Simulation of a space environment requires that gravity 
and friction be negligible. A large granite table serves as 
the foundation for the laboratory set-up because precision 
flatness and smoothness can be achieved. The mainbody and 
flexible arm float on air pads that skim over the surface of 
the granite table. 

The ta^^le is 6' X 8' X 10.5" and rests on a castor stand. 
There are set-screws on the metal stand to achieve a level 
platform. The surface is a Laboratory Grade A .001" finish. 
The stand has wheels to allow mobility. The table weighs 
7550 lbs. 
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B. MAINBODY 


The overall size of the mainbody is constrained by the 
dimensions of the granite table. In order to have adequate 
room to maneuver the system when slewing, a 30 inch diameter 
was chosen. 

The base of the mainbody is a simple 7/8" aluminum disk 
with a 15 inch radius (Figure 2). The mainbody was machined 
at Naval Postgraduate School by the Aeronautical and 
Astronautical Engineering department. The top side of the 
mainbody has seventy two 1/4 inch #20 drilled holes 
positioned in a radial pattern to allow for a wide 
assortment of attachments as the experiment grows. 
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To achieve a realistic simulation of actual spacecraft 
conditions, the moment of inertia of the mainbody is assumed 
to be roughly three times the moment of inertia of the 
flexible arm. From this specification, the moment of 
inertia of the mainbody can be calculated: 


^mainbody ^ ^ ^arm assembly'^ 

^mainbody = (3) X 25,000 Ibw-inches^ (1) 

^mair^ody = 75,000 IbW-jncheS^ 


The total mass for the mainbody can be calculated from 
its moment of inertia: 


■^mainbody ~ ~ 2 ^mainbody^inainbody 
^mainbody = ^67 IblTl 


Where r„,. = 15 inches 

ma 1 nbody 


( 2 ) 


For practicality in handling, the mass of the entire 
mainbody should be under 500 lb. This means that the 
inertia ratio of flexible assembly/mainbody will be more on 
the order of 2 to 1, or an of approximately 50,000 lb- 
in^. This gives a total moment of inertia for the entire 
system of approximately 75,000 Ib-in^. 

The momentum wheel and motor assembly will add 
approximately 30 lb to the system. From the above 
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calculation, it is obvious that additional weights must be 
added to the base to achieve the desired moment of inertia. 

The mass of the base of the mainbody can be calculated as 
follows: 


M, 

M, 


base 


base 


P ^ ^base^base 
= 59.376 Ihm 


(3) 


wherep = .096 Ib-in^ 
^ base = inches 

h base = inches 


From the above, the moment of inertia of the base of the 
mainbody can be calculated by: 

1 

^mb ~ ^mb^mb (4) 

^mb ~ 6,67 9.80 Ibm-inches'^ 


C. FLEXIBLE ARM 

The flexible arm (Figure 3) consists of two assemblies of 
aluminum and steel construction bolted together rigidly in a 
right angle elbow. Mass intensifiers are connected to the 
thin (.16 inch) aluminum bar to increase the moment of 
inertia of the arm assembly without significantly increasing 
the stiffness. Each mass intensifier (Figure 4) consists of 
two parts which are bolted together on each side of the 
aluminum. The tapered design of the mass intensifier 
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reduces the contact between the steel and the aluminum. 
Minimal contact will ensure greater flexibility. 
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Figure 4 Mass Intensifiers. 


The elbow joint consists of two aluminum right angle 
braces connected to a circular aluminum disk. The end joint 
has an aluminum wedge bolted to a circular aluminum disk. 
This wedge shape allows free motion of the endpoint. The 
end of the flexible arm will house a reflector and will be 
the location for end point sensing. 

The total mass of the flexible arm assembly can be 
calculated as follows: 

mass of aluminum rod^ = .411 Ibm 
mass of aluminum rodg = .380 Ibm 
mass of each intensifier pair = 1.056 Ibm 

Total mass of arm^^ = 9.66 Ibm 
Total mass of arm^ = 8.59 Ibm 

mass of elbow joint = .244 Ibm 
mass of end joint = .183 Ibm 

Total mass of flexible assembly = 18.68 Ibm 
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The center of mass of the flexible arm assembly is: 


Center of mass of arm. 


= 1 : 




CM. 


arm* 


^arm* 

12.96 inches 1 


Center of mass of armg = 


r 


E 


Af. 


CM. 


arnin 


= 11.66 Inches j 


Center of mass of arm assembly ' 

assembly 
^y) assembly ~ 


^crrr^ cm 

M 


CM^r^ ( X) + M^^^ CM^rm (^) 

arm% ariOt ' ' a.Tmg ATWg ' 


Warm.O^arm. '^V) (y> 


rm* '*’^a rm* 


C^(>^)aaa«mt.iv- = 19.29 incheS 
CMf.y) assembly = 5.55 inches 


(6) 


The moment of inertia of each arm can be calculated about 
the center of mass of each arm. Each mass intensifier and 
the elbow and end joint braces are computed as point masses. 
The moment of inertia for the flexible arm assembly can be 
calculated about the origin O by the parallel axis theorem. 
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) = 416.10 Ibm-inches^ 

(J™, ) = 324.03 Ibirt-inches" 




M. 


where d - (15+12.96) = 27.96 inches 


(J,^_ ) = 7,967.9 2 Ibm-inches^ 


(7) 


() 


^ ^armg 


where d = v^(15 + 27.96 ) ^+(11.66 ) ^ = 42.85 inches 


(J... ) = 16,096.32 Ibm-inches^ 
*0 


■^ar/n assemhiy^ ^ 

assemliy^ = 24,064.24 Ibm-inches^ 


D. MOMENTUM WHEEL 

A motor driven momentum wheel will be used to apply a 
torque (the change in the angular momentum) for slewing the 
mainbody. 

For preliminary calculations, a one foot diameter 
aluminum disk with one inch thickness can be used for the 
dimensions of the momentum wheel. 
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The mass of the momentum wheel can be calculated as: 



P ^ ^mw 

Where r^^ 

=6.00 inches 

^mw 

=1.0 inches 

P 

. .096 


inches 


= 10.857 Ibm 


(8) 


The moment of inertia for the momentum wheel is: 


I„ 


where 


I. 

I. 



6.00 inches 
10.857 Ibw 
195.43 Ibw-inches^ 
8.094 oz-inches-sec^ 


(9) 


E. TORQUE MOTOR 

The momentum wheel must be driven by a motor that will 
deliver the required change in angular momentum to the 
system. A specification for the system design, based upon 
existing space platforms, is that the mainbody accelerate at 
the rate of 12 degrees per second^. 
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From this specification, the torque of the motor can be 


determined to be: 


^ syscem^ Eyscem 


(3000oz-irches-sec^) (.209 


rad 

sec^ 

a. 


] = (8.094oz-inches-sec^) (a ) 


= 80 


Torque of momentum wheel 


rad 

sec^ 


( 10 ) 


= 6 50 oz-inch 


From the above calculation, the momentum wheel must apply 
a 650 oz-inch (4.59 N-m) torque to provide the desired slew 
rate for the mainbody. 

The motor selected was a DC servo-disc motor. In the 
servo-disc motor, the ironless armature is constructed from 
several layers of copper conductors in a flat-disc 
configuration. This allows for low inertia and fast 
acceleration. Additionally, the smaller size can be easily 
mounted on the mainbody. Table I outlines the motor 
characteristic. 
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TABLE I PMI SERVO-DISC MOTOR 

Motor Performance 

Units 

TR16M4C 

Peak Torque 

oz-in 

5307.2 

Continuous Stall Torque 

oz-in 

498.4 

Peak Current 

Amps 

100.8 

Peak Acceleration w/o Load 

KRADS/sec^ 

63.2 

Horsepower 

HP 

1.4 

Torque 

oz-in 


Speed 

RPM 

3000 

Power Output 

Watts 

1049.5 

Terminal Voltage 

Volts 

128.7 

Current 

Amps 

9.55 

Torque Constant (KT) 

oz-in/ 

52.77 

Motor Weight 

Lb 

17.50 

Moment of Inertia 

oz-in-s^ 

. 084 

Cost (includes tachometer) 

ea 

$1229.0 


The motor will be aligned in a vertical position on top 
of the momentum wheel to deliver a torque about the vertical 
(pitch) axis. The motor housing includes an attachment 
shaft for fitting with the momentum wheel. 
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F. SENSORS 


Initially, the angular position and the angular rate of 
the mainbody and the speed of the momentum wheel will be 
sensed. At later time, the angular position of the 
reflector with respect to the endpoint, and the position of 
the endpoint with respect to the mainbody, will be measui d. 

The angular rate of the torquer mooor and momentum wheel 
will be sensed by a tachometer which is attached to the 
motor. Table II outlines the tachometer characteristics. 


TABLE II PMI TACHOMETER 

Tachometer Performance 

Units 

JR16M4CH 

Maximum Speed 

RPM 

4000 

Bidirectional Tolerance 

%V 

±1 • 5 

Tachometer Weight 

lb 

2.2 

Moment of Inertia 

oz- 

. 005 


in-s^ 



The angular rate of the mainbody will be monitored by a 
rate sensing gyroscope. The gyroscope produces an analog 
output voltage proportional to the angular velocity about 
the sensing axis ''in this case, the Z axis at the origin). 
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When a rotation exists, the Coriolis forces transfer 
momentum to the perpendicular plane and cause a bending of 
the solid state sensing elements. These elements will 
vibrate 180 degrees out of phase with angular motion but in 
phase when linear acceleration or vibration occurs. The 
integrated electronics produce an angular motion response 
output with no linear motion incorporated. A single axis 
gyroscope with a range of +30 degrees per second was 
selected. Table III shows the rate gyroscope 
characteristics. 


TABLE III WATSON RATE GYROSCOPE 

Rate Gyroscope 

Units 

ARS-C121-1A 

Output 

VDC 

0 at zero angular rate 
+10 at full scale 
angular rate 

Power Supply 

VDC 

+ 15 + 5% 20 mA 
maximum 

Range 

degree 

/sec 

+ 30 

Weight 

Oz 

4 

Cost 

ea 

$ 802.00 


Angular position for both the mainbody and the endpoint 
reflector will be measured by rotary variable differential 
transformers (RVDTs). RVDTs measure the angular 
displacement of rotating elements by producing a voltage 
whose magnitude varies linearly with the angular position of 
the shaft. The RVDT is constructed on precision ball 


























bearings to minimize friction and noise. The best operating 
range of a RVDT is + 40 degrees, which is well within the 
confines of this experiment. The characteristics of the 
RVDTs purchased are in Table IV. 


TABLE IV SCHAEVITZ ROTARY 
VARIABLE DIFFERENTIAL TRANSFORMER 


RVDTs 

Units 

R30D 

Linearity 

% of 

+ 3 0° : .25 


range 

+ 40° : .50 



+ 60° : 1.5 

Moment of Inertia 

Oz-in- 

8.5 X 10’'^ 


sec^ 


Maximum Torque 

Inch-oz 

. 019 

Maximum Load 

Lb 

Radial = 8 



Axial = 10 

Weight 

Oz 

1.87 

Cost 

ea 

$438.00 



Endpoint sensing has yet to be determined. The 
accuracies required (+ .1 mm) over the large span of control 
(+ 5 cm) make most conventional endpoint sensing techniques 
inadequate. CCD cameras would provide the required coverage 
area, however they would give only 1 mm accuracies and the 
































signal processing equipment required to integrate CCD 
cameras into the control loop would be expansive. Greater 
accuracies can easily be achieved by off-the-shelf 
capacitance sensors, however the coverage area would be much 
smaller than that required by this experiment. 

A promising devise for lateral position sensing is the 
electro-optical photodetector. Using a fixed-position laser 
or LED light source, a mirror system and a photodiode 
receiver collect the reflected light and provide an output 
signal proportional to the position. 

Standa^-d off-the-shelf units can be purchased through 
several sources. The unit would consist of an analog 
optical position monitor (approximate cost $3,500) and a 
detector/cable assembly (approximate cost $400) . 

G. AIR PADS 

The mainbody and the flexible arm will be supported by 
five aluminum air pads, each capable of supporting a 60 lb 
load. Each pad is fed by 80 PSI of air, which is regulated 
by two regulators. One regulator feeds the three mainbody 
air pads, and one regulator feeds the elbow and endpoint air 
pads. 

On the mainbody, three 1.5 inch diameter mounting holes 
are located 120 degrees apart three inches from the outer 
edge. The air feed is on the side of each pad. A top brace 
will be fashioned and mounted to the mainbody to hold the 
air pads in place. Three air pads for the mainbody have 
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been chosen as the optimum number for controlling the 
mass. An odd number of pads decreases the likelihood of 
leverage effects due to uneven load distributions. A 150 
psi regulator will control all three mainbody air pads. The 
150 psi regulator was chosen to allow upgrades to larger air 
pads at a later date. 

H. AIR BEARING 

The entire system will be constrained to rotational 
motion only by an air bearing mounted to the center of the 
mainbrtdy from above. Translational motion will be 
constrained in the initial phase of the experiment so that 
pitch effects can be isolated. At a later time, the air 
bearing can be removed so that robotics experiments can be 
performed. The air bearing (Figure 5) consists of an H-type 
bearing which is capable of carrying both a radial and an 
axial load. As the mass of the mainbody is increased, the 
air bearing can be used to support some of the load. A 
mounting bracket will be attached to the side of the granite 
table and will be connected to the mainbody from above. 
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Figure 5 Air Bearing. 


I. COMPUTER SYSTEM 

The AC-100 computer system, manufactured by Integrated 
Systems, Inc., has been purchased to provide data 
acquisition and results analysis for the experimental set¬ 
up. The AC-100 is an automatically-programmable real-time 
control system that will allow graphical specification and 
simulation of real-time systems. This will allow the 
experimental set-up to be fully monitored real-time from the 
workstation. 

The AC-100 system consists of three parts; the 
workstation, the AC-100 software, and the AC-100 controller. 
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The workstation includes a DEC 3100 computer, VMS operating 
system, 8 MB RAM memory, two 105 MB hard disks, a 19 inch 
color monitor and an Ethernet interface. PWS003 modeling, 
simulation and implementation software are combined with a 
VAX C and VAX Fortran compiler. The controller has 16 
channel inputs and 10 channel outputs in addition to 32 
parallel digital inputs/outputs, eight encoder inputs, and 
one event trigger input. 
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III. EXPERIMENTAL ANALYSIS 


A. MODAL ANALYSIS 

Experimental modal testing can be conducted to determine 
the nature of the vibration response and to verify the 
analytical models. The Hewlett-Packard VISTA Data 
Acquisition Package was used to collect information 
concerning the time and frequency response of the flexible 
arm assembly in a free-free mode. 

The frequency response function method of modal testing 
requires that the input excitation and the output response 
be measured simultanously to determine the system function 
(Figure 6). 



The VISTA program utilizes a dynamic signal analyzer, 
which is a Fourier transform-based instrument, to process 
the measurements of the structures frequency response 
(Hewlett- Packard, 1986, pp. 11). 


27 










The input used was a single-point excitation impact 
hammer. The hammer is the most simple way to excite a 
structure into vibration. The exciter consists of an 
impactor, with several different tips and heads which will 
allow some flexibility in setting the force level ranges for 
testing different objects (Ewins, 1984, pp. 102-104). The 
hammer contains an ICP quartz force sensor mounted on the 
striking end of the hammer (PCB Piezotronics, 1983, pp. 1). 
The impact force is transferred via the sensing element into 
an electrical signal which can be evaluated. 

The hammer impulse consists of a nearly constant force 
over a wide frequency range and therefore can excite all 
resonances in that range. The size of the hammer, along 
with the hammer tip material and velocity, determine the 
amplitude of the force impulse. The frequency content of 
the energy applied to the arm is a function of the stiffness 
of the contacting surfaces and the mass of the hammer. The 
hammer tip selected affects the force impulse, and therefore 
the frequency response. A hard, steel tip was selected to 
deliver a short pulse duration and higher frequency content. 

The flexible arm assembly was supported from above by 
small soft elastic cords in an attempt to simulate the free- 
free mode (Figure 7). An ICP accelerometer was positioned 
in the XY plane at the end of each arm. When the 
accelerometer vibrates, an internal mass applies a force to 
the crystal element which is proportional to the 
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acceleration. Using Newton's Law, the resulting force 
output can be measured. 



Figure 7 Flexible Arm Assembly Support. 


The assembly was struck with the hammer in various places 
throughout the assembly. The resulting frequency response 
showed consistent modes of varying amplitudes, depending 
upon where the structure was struck. 

Figure 8 shows the time and frequency spectrum for the 
impact hammer as struck about midway on arm A. Figure 9 and 
Figure 10 show the time and frequency measurements for each 
arm as a result of the impact. 
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Figure 9 Arm A Time and Frequency Response. 



Figure 10 Arm B Time and Frequency Measurement. 
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The finite element analysis was evaluated for the free- 
free condition, and the mode frequencies were computed. 

When compared with the experimental results, a fairly good 
correlation was found between the two, with the exception of 
the first mode. Table V shows the difference between the 
frequencies obtained from the GIFTS program and the 
experimental values. 



TABLE V FINITE ELEMENT ANALYSIS 

VS. EXPERIMENTAL RESULTS 

MODE 

GIFTS 

Hz 

EXP.AVG 

Hz 

DIFFER¬ 

ENCE 

Hz 

1 

. 548 

1.125 

-.577 

2 

2.087 

2.249 

-.162 

3 

3.057 

3.375 

-.318 

4 

6.483 

6.093 

. 390 

5 

8.133 

8.312 

-.179 

6 

13.424 

12.313 

1.111 

7 

15.667 

15.125 

. 542 
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IV. THEORETICAL ANALYSIS 


A. EQUATIONS OF MOTION 

The system equations of motion can be derived as a 
mathematical model for the exoerimental set-up. This model 
is based upon given physical laws which govern the system 
behavior, such as force-deformation relationships and the 
laws of motion. The equations of motion wil? identify the 
mass and stiffness system parameters. 

In a "lumped-parameter” system, the equations are a 
function of time alone, which are ordinary differential 
equations and are easy to solve (Meirovitch, 1990, pp. 45- 
53). In the "distributed-system-parau.eter" system, the 
equations are functions of both spatial coordinates and 
time, and are therefore partial differential equatic .s which 
can not be easily solved. Therefore, a discret model must 
be devised to model the distributed system. 

A flexible system is a non-linear system. In other 
words, the response of the experimental set-up to different 
excitations can not be obtained separately and combined 
linearly. Nonlinearity may present complications in 
modelling and control. A careful ccmparison of the 
analytical verses experimental results will identify 
nonlinear variations. 
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The equations of motion for a flexible structure can be 
written as a function of spatial variables, which can be 
assumed cjs known, and unknown time variables. This will 
result in a discrete model containing ordinary differential 
equations. The solution can be represented as spatial 
eigenfunctions multiplied by time-dependant generalized 
coordinates. 

Lagrange's equation facilitates writing the equations of 
motion for a flexible multiple degree of freedom system. 
Additionally, Lagrange's technique produces an entirely 
scalar process using the quantities of potential energy, 
kinetic energy, and work in terms of a generalized 
coordinate system (Meirovitch, 1967, pp.47-59). 

The Lagrange Operator is: 


d , dL 

d c ' 



(It) 


where L is the Lagrangian: T - V 
T = System kinetic energy 
V = System potential energy 
q- = System generalized coordinates 
(b = SyF.:e . generalized force 
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By using this method of analysis, the complicated dynamic 
problem can be reduced to a scalar integral that does not 
depend on the coordinates used. 

The kinetic energy of a distributed system is: 

T=f^p(V^)^ ds (12) 

s 


where p is the material density, is the velocity of a 
generic point of the system in inertial frame 


The origin O for the experimental set-up was chosen to be 
the center of the mainbody. This is the point where the 
mainbody is attached to the air bearing (Figure 11). This 
origin O is therefore fixed in the reference frame of the 
table. 

The value of can be determined from the following 
formula: 


where o) = 6 ^ 


(13) 
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The position vector from the origin O to any point on the 
arm can be described by: 

?o-a = X i + j 
^o-b = (R+Li+W,) i + + y) j 

where = length of Arm A 
R = radius of the mainbody 
Wgg = the vertical endpoint displacement of Arm A 

thus 

to X fo .3 = X [ iXi+W^j) ] = -Wjbi+X^j 

00 X = 6ic X [(R+Li + f^) 1 + (y+Wag) j] (14) 

= 1 + (i?+Li+f4)6 j 

Using the function Wg(x,y,t) to represent the position 
of any arbitrary point on beam A, and to represent 

the position of any arbitrary point on beam B: 

W, =^\f^,{x)g^{t)*^^,ix)g,{t) 

= ^2 {L^} g2 i t) 

=Pi g^ ( t) +P2 g2 ( t) 

where t is a function of spatial variables only and 
q(t) is the time dependant general coordinate 
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The subscripts 1 and 2 refer to the first and second 
inodes of the structure. Only the first two inodes have been 
incorporated in the spatial equations, but others could 
easily be added at a later time. The displacement in X 
direction for Arm A has been ignored as insignificant 
compared to the other displacements. 
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Substituting this new variable into the equation, the 


kinetic energy of the system can be calculated to be: 


KEsystein= — I^ 


R*L, 


^0 R ^0 

R*Li L 2 

[pJpiPjdy+p I 4r^_j4r^2<*<^+p|4r^jilfi,2dy] 

OR 0 

L 2 R*L^ L 2 

[-|/p2c?y+^ / 4fLd>^+-^/4fLc?y] 

0 R ^ 0 

R*Li L 2 

R ^ 0 

R*L^ L 2 

[PP 1 P 2 L 2 + P I taj4'a2^+p/’ifijllffc2C?y] 

R 0 

L 2 

[p (i?+Li)L2pi + P / X}\f^jdx-pjY}^j^^dy] 

R 0 

K+Li L 2 

[pp2L2+-| J tLdx+^JtLciy] 

i? ^0 

i? * 

[p (i?i+L^)L2P2 + P I xilr^^dx-pjyilf^^dy] 

K 0 

Lj 

[-^ j x^dx+-^ Jy^+(i?+Li) ^dy] 

[pj (i?+Lj i|;^^dy] q^B^+[pj (i?+Li) 

0 0 

J^2 ^ 

tp/p2'l'fci^y] ^1^72^ " fp/p2tfc2C?y] 

0 0 

[p|p2rdy] g2^^ + [p/Pii"c?y] 

0 0 
L 2 1*2 

ip/p,* t2<iy] gi^zP - [p/PitwCfyl 


e^+ 

gfe^- 

<7iQ-2p" + 

<7iP + 

<72p + 

g20^ - 
gzcjjP + 

g^e^ - 

<7igi0 


( 16 ) 
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The potential energy of the system is the flexural strain 
energy. To calculate the potential energy of each beam, the 
following equation can be used: 


R*L^ 


y-l / EI( 


dx^ 


) ^dx + 


1 

2 



dX 

dy^ 


) ^dy 


( 12 ) 


Substituting the expressions for W^, and into this 

equation, the total potential energy of the system is 
calculated. 

The rotary inertia effect and the axial deformation 
effect on the potential energy have been ignored. The cross- 
sectional dimensions are small compared to the length of the 
bar. 
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The Lagrangian can now be written in the following form 


( 18 ) 




R ^0 

R*L^ 1/2 

[■|/P?c?y+-| f ^hdx-t^-^fiilidy] qffi^ + 

0 J? 0 

Z/^ R*Z/y Z /2 

[pJpjPjdy+p f ^ 

0 j; 0 

£•2 R*Z/y £>2 

[-|/P2C?y+-| / '¥l2dx*-^j^l2dy] q|6^ + 

0 ^ 0 

^♦Lj L 2 

+ f 'lf^idx*-^J‘}lfijdy] c*f^ + 

1/2 

CPPlP2-t'2+P / 'l'flJ»l»a2CiX+p/ti,2'l'i,2'^y] + 

R 0 

R*L. 1/2 

[p(/?+Lj)L2Pi + p J Xilr^jdx-pJyilr^jdy] dj6 + 

J? 0 

R*l^ L 3 

R*L>, Li 

[p (i? + Li)L2P2'*-p f + 

R 0 

[p|(i?+Li)i|rj^jdy] gi6^+[pj (iJ+Lj il/^^dy] - 

0 0 

Li Lj 

tp/+wP2C(K] '^^i< 72P - p/P2'l'2,2C?y] <724^:6 + 

0 0 

[p/P2i'cfy] <72®^ + [p/Pi5^c?y] <7iP^ - 

0 0 
Li Lj 

tp/H'i,2PiC?y] Qi<^2^ - [p/ikwPiC?yJ gi^iP - 
0 0 

R*L. Li 


R*L^ 


[£•7 f ( 




) dx+Elj" { 


d^i|r^j d^il;^ 


R y//^ -Q dy^ dy 


- “ '^-)dy] g^g^ 


dx^ dx^ 

ff*L.^ 5 ^ 
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Performing the Lagrangian Operator, the resulting matrix 
is as follows: 

Mil M12 M13 Kll K12 K13 Qi 0 

M21 M22 M23 gs + K21 K22 K23 ^2 + ^2 = 0 

M31 M32 M33 0 K31 K32 K33 0 

where : 


Mil 

Ml 2 

Ml 3 

M21 

M2 2 

M2 3 

M31 

M3 2 

M3 3 



( 20 ) 
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i { dy^ 

ElT^^dx^Ell^ ^ 

•|j dx^ dx^ J dy^ dy^ 

^dx.Ejf-£!iM ^ dy 

•<_ d>r^ d>r^ J dy^ dy^ 


Elf () ^dx^EI f ( —^ ) 2dy 

i i dy^ 

K23 = K31 = ii:32 = i(:33 = 0 


Lj 

Vi = -lpfpldy->-p I ^iicbc-t-p^i^l^dy] - 

0 R 0 

R+I.. Lj 

[p/PiPjCfy+P / - 

0 ^ 0 

1-2 j^2 

[pj (i?+Li) ffcjdy-[p^Piydy] 6^ - 

0 0 

tp/p2’l»i.iCfy-p/piT|>j,^dy] - 

0 0 

fp/p2ti,jc?y] gjP + 

0 

tp/Pi4’w<^y] 

0 

t>2 i?*Lj ^2 

l^^-tp/PiPsCfy+P / ’I'aitarCLx'-^p/twth^rfy] diP" - 

OR 0 

1.2 Lj 

[pfpidy*p J ^)»Ldx+p|i|»|2dy] g^P^ - 

0 ft 0 

[pJP2ydy+pJ(i?+I,i) i|r^^dy] P^ - 

0 0 

iZpfp^'I'bsdyl <^2® - 
0 

fp/Pi’l'w<^y"^p/'l'wp2c?y] - 

0 0 

ip/Pi+wcfyi <7iP - 
0 

tp/P2t/>2C?y^ ^2® 
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—i —1 < 

V3 = [2pj'pldy+2p J i|r^idx+2pji|fwdy] + 

0 A 0 

R*L 2 

[2pj'pjdy+2p J \(r^jdx-+2pj‘i|»|,;,dy] ^1(^16 + 

0 J? 0 

ZpJ ^*X*J i»2 

[2pj'PiP2dy+2p J i|r^i4'a2tix+2pji|i^_,i(r^^dy] + 

0 0 

[2pJ(i?+I.i)i|j^jdy+2pj‘Piydy] 4ri6 + 

0 0 

[ 2 pJ(i?+Lj 4r^2dy+p|P2ydy] <^26 + 

0 0 

Z»2 ^♦X.j Z»2 

[2p|piP2dy>2p / 4'a2’lf«2dx-+2p|iltijT|ri,^dy] gi(3r2e 4 
0 ^ 0 

R*L^ 

[p/Pirfy+p / tLci>f+p/4'Lc?y] + 

0 /e 0 

R*Li-i Lt 2 

[2p|p,p2dy+2p / <l»a2tai<dy+2p/^liiji|»fc^dy] g^gjP f 
0 0 0 

X<2 ^♦Xn X>2 

[pJPjdy+p J i|»Lcix+pJ+lidy] gffi + 

0 * > 

I, X, 

[2p|(i?+Li)»|>2,jdy+2pJP2ydy] gjP + 

0 0 

1-2 ^2 ^ 

[2p|(i?+I,i)\lfjb2C?y+2p|P2ydy] g^e - [pJPjilfwCfy] (Jigz - 

0 0 0 

L, L 2 X, 

[p/Patbirfy-p/Pit^dy] (jr^gj -[pJPzt^^dy] - 
0 0 0 

[p/P2’l'w‘^y] '72<3^2 ■^tp/Pi’l'w‘*'5 <71 <^1 " 

0 0 
X-2 1*3 

[p/Pi’i'w'^y] - [p/+t.jPic?yi 
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The spatial variables can be determined in several ways. 


The assumed-method (Kane, 1983, pp. 318-319) makes the 
approximation that the deflection of each flexible link can 
be expressed as a sum of a finite number of spatial mode 
shapes that satisfy the given boundary conditions. For a 
uniform cantilever free-free or fixed-free beam, the 
components of the spatial variables have been calculated by 
numerous sources. 

The experimental set-up, however, is not a simple 
cantilever/free-free continuous beam construction and 
therefore would not conform to the above simplification. 
Therefore, a finite element analysis of the flexible arm was 
conducted to determine the first two mode shapes. The 
finite element analysis program GIFTS was used to create a 
model of a continuous aluminum bar. Point mass loads were 
introduced at the position of the mass intensifiers and the 
aluminum connecting elbow and end braces. The structure was 
divided into forty nodes. 

For a fixed-free representation, the position 
was supported in the finite element analysis model. From 
the GIFTS analysis, unit deflections for each mode shape 
could be determined. Only the first two modes were 
analyzed. Each unit was arbitrarily determined to be one 
inch. By adding the unit deflection to the original 
position, a polynomial could be fitted for each mode shape 
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(Figui.es 12 & 13) . Therefore, the spatial variables are as 
follows: 



= - .0057 

- 

.0345 

X + 

. 0679 

'I'm 

= .0009 

y2 

.3233 

Y + 

27.72 

'I'a2 

= -.0025 

+ 

.2462 

X - 

.4822 

'J'm 

= .0074 

y2 ^ 

.1958 

Y + 

25.17 


The vertical deflection of Arm A endpoint can be 
determined by evaluating this function to be: 


= -4.77 
- 4.26 



RRI I B 

psini 

RRM 1 

MODE ] 
RRM B 

MODE ] 
RRM R 


Figure 12 MODE 1 FINITE ELEMENT ANALYSIS 

VS. POLYNOMIAL FIT. 


(24) 










Figure 13 MODE 2 FINITE ELEMENT ANALYSIS 

VS. POLYNOMIAL FIT. 


Performing the assigned integration results in the 
3olution of the matrix variables as foliowu: 


Mil = 202.84 
il2 = -44.65 
M13 =-316.25 
M21 = -44.65 

M22 = 179.05 (25) 

M23 = 170.81 

M31 =-316.25 

M32 = 170.81 

M33 = + 168.14 
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= - 153.18 -18.57 q^^' - 

398.68 6 ^ + 39.84 q^^ - 
19.72 ^20 - 22.08 

= -155.7 8 g^^^ - 52.50 g 26 ^ - 

363.07 6 ^ - 35.94 (^26 + 

.4 <^^6 + 20.12 gi 6 - 
17.97 g20 

V 3 = 305.00 gj^j^ + 306.36 (26) 

+ 13.50 + 13.50 + 

71.0.32 <^36 +711.11 ^ 2 ^ + 

153.18 gi0 + u.3.50 g2g20 + 

152.50 g^0 + 793.90 g^S 
711.11 g20 - 19.72 g 3 g 2 - 
39.84 q^q^ - 11.91 q^ - 
17.97 g^g, - 22.08 g^g^ + 

22.08 q( + 22 .OBqj^q, 


Kll = .71 

K12 = .435 

K21 = .435 

K22 = 1.25 


(27) 
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B. NATURAL FREQUENCY 


When placed in motion, the oscillation of a system that 
takes place will be at the natural frequency u^, which is a 
property of the system. Damping in small amounts will have 
little effect on the natural frequency and has been 
neglected in these preliminary calculations. For a single 
degree of freedom linear system, the natural frequency can 
be determined from the equations of motion by: 


where K is the stiffness matrix 
M is the mass matrix 
For the experimental set-up; 


E = 10 X 10* Ibm/in"^ 

I = 2.0345 X 10'^ in^ 

= .381 Ibm/in 

Ignoring the nonlinear terms, solving the equation for 


= . 0460 Hz 
(*>2 = • 1004 Hz 
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V. 


CONCLUSIONS AND RECOMMENDATIONS 


Based on this work, the following conclusions are 
reached. 

A. SYSTEM DESIGN 

The experimental set-up for Phase I study of spacecraft 
control/structural interaction has been designed. All the 
actuators and sensors have been selected except the end¬ 
point displacement sensing of the arm. The mainbody and the 
flexible arm have been fabricated. 

B. THEORETICAL ANALYSIS 

The equations of motion for the experimental model have 
been derived and natural frequencies determined. 

C. EXPERIMENTAL ANALYSIS 

The natural frequencies of the flexible arm has been 
determined experimentally and compared with analytical 
predictions obtained by using the GIFTS finite element 
analysis program. The experimental and analytical results 
are in good agreement except the first mode. 

D. RECOMMENDATIONS 

1. The analytical model of the experimental set-up needs 
to be improved by full representations of flexible modes, 
including control laws, and considering nonlinear effects. 
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2. Future modal analysis should incorporate a filter to 
reduce the effects of noise. Force windows and exponential 
windows are available in the data acquisition system for 
such a purpose. 

3. Modal analysis using the VISTA data acquisition 
package and IDEAS modal analysis software should be 
completed on the structure once it is connected together and 
mounted on the granite table on airpads. Natural 
frequencies and mode shapes should be then determined. 
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